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1.  INTRODUCTION 

The  South  China  Sea  (SCS),  Yellow/East 
China  Sea  (YES),  and  Japan/East  Sea  (JES)  are 
major  East  Asian  marginal  seas  (EAMS).  The 
complex  topography  includes  the  broad  shallows 
of  the  Sunda  Shelf  in  the  south/southwest  of  SCS; 
the  continental  shelf  of  the  Asian  landmass  in  the 
north,  extending  from  the  Gulf  of  Tonkin  to  the 
YES;  a  deep,  elliptical  shaped  SCS  and  JES 
basins,  and  numerous  reef  islands  and 
underwater  plateaus  scattered  throughout  the 
basin  (Fig.  la).  The  shelf  that  extends  from  the 
Gulf  of  Tonkin  to  the  YES  is  consistently  near  70 
m  deep,  and  averages  150  km  in  width. 

The  EAMS  is  subjected  to  a  seasonal 
monsoon  system.  From  April  to  August,  the 
weaker  southwesterly  summer  monsoon  winds 
result  in  a  wind  stress  of  just  over  0.1  N/m2.  From 
November  to  March,  the  stronger  northeasterly 
winter  monsoon  winds  correspond  to  a  maximum 
wind  stress  of  nearly  0.3  N/m.  Recent 
observational  studies  show  that  the  EAMS  is 
energetic  and  has  multi-eddy  structure.  For 
example,  the  SCS  synoptic  eddy  structure  was 
identified  in  May  1995  using  the  airborne 
expendable  bathythermograph  (AXBT)  data  (Chu 
et  al.,  1998a),  the  eddy  spatial  and  temporal 
scales  in  the  YES  were  identified  using  the  Navy's 
Master  Oceanographic  Observational  Data  Set 
(MOODS)  during  1929-1991  (Chuetal., 

1997a,b),  and  the  seasonal  JES  multi-eddy 
structure  from  a  composite  analysis  on  the  U.S. 
National  Centers  for  Environmental  Prediction 
(NCEP)  monthly  SST  fields  during  1981-1994 
(Chu  et  al.  1998b). 

The  fundamental  scientific  issues  are:  What 
are  the  dynamical  balances  controlling  mesoscale 
eddy  variations  in  EAMS?  What  are  effects  of 
surface  wind  and  thermohaline  forcing,  coastline 
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geometry  and  topographic  slope  on  the  coastal 
dynamics  including  wave  processes,  current 
instabilities,?  What  is  the  role  of  coastal  air-ocean 
coupling  on  oceanic  and  atmospheric  dynamics  of 
mesoscale  eddies? 

To  study  these  problems,  a  coastal 
atmosphere-ocean  coupled  system  (CAOCS)  was 
developed  at  the  Naval  Postgraduate  School.  The 
model  domain  (67°-142°E,  7°S-55°N)  covers  the 
whole  EAMS  and  surrounding  land  and  islands. 
The  surface  fluxes  of  water,  heat  (excluding  solar 
radiation),  and  momentum  are  applied 
synchronously  with  opposite  signs  in  the 
atmosphere  and  ocean.  Flux  adjustments  are  not 
used. 

2.  CAOCS  FOR  EAMS  NOWCAST/FORECAST 

2.1  Ocean  Component 

2.1.1.  Numerics 

The  CAOCS  ocean  component  (POM)  is 
a  time  dependent,  primitive  equation  circulation 
model  on  a  three  dimensional  grid  that  includes 
realistic  topography  and  a  free  surface.  From  a 
series  of  numerical  experiments  (Chu  et  al., 

1998a;  1999a,b;  2000b),  the  qualitative  and 
quantitative  effects  of  non-linearity,  wind  forcing 
and  lateral  boundary  transport  on  the  EAMS  are 
analyzed,  yielding  considerable  insight  into  the 
external  factors  affecting  the  regional 
oceanography.  We  use  a  rectilinear  grid  with 
horizontal  resolution  of  0.25°  and  23  non-uniform 
vertical  a  levels.  The  model  uses  realistic 
bathymetry  data  from  the  Naval  Oceanographic 
Office  Digital  Bathymetry  Data  Base  with  5 
minutes  resolution  (DBDB5). 

2.1.2.  Open  Boundary  Conditions 

Closed  lateral  boundaries,  i.e. ,  the 
modeled  ocean  bordered  by  land,  were 


defined  using  a  free  slip  condition  for  velocity  and 
a  zero  gradient  condition  for  temperature  and 
salinity.  No  advective  or  diffusive  heat,  salt 
or  velocity  fluxes  occur  through  these  boundaries. 

At  open  boundaries,  we  use  monthly 
varying  volume  transports  from  a  global  inverse 
model  (Chu  and  Fan,  2000)  for  the  open 
boundaries  (Fig.  la).  When  the  water  flows  into 
the  model  domain,  temperature  and  salinity  at  the 
open  boundary  are  likewise  prescribed  from  the 
climatological  data.  When  water  flows  out  of  the 
domain,  the  radiation  condition  was  applied, 

^-(a,s)+un^-(a,s)  =  o 

at  an 

where  the  subscript  n  is  the  direction  normal  to 
the  boundary. 

2.1.3.  Initial  Conditions  and  Model  Initialization 

Before  coupling  to  the  atmospheric  model 
that  is  MM5,  the  POM  model  was  integrated  for 
four  years  and  four  months  from  zero  velocity, 
and  January  climatological  temperature  and 
salinity  fields  and  forced  by  monthly  mean 
surface  wind  stress  from  the  Comprehensive 
Ocean  and  Atmosphere  Data  Set  (COADS)  and 
by  the  restoring-type  surface  heat  and  salinity 
fluxes  which  are  relaxed  to  the  surface  monthly 
values.  The  final  states  of  temperature, 
salinity,  and  velocity  were  taken  as  the  initial 
ocean  conditions  for  1  April  1998. 

2.2.  Atmospheric  Component 

2.2.1.  Model  Description 

We  take  Version-3.4  of  the  Fifth- 
Generation  Penn  State  University/NCAR 
Mesoscale  Model  (MM5)  as  the  atmospheric 
component  for  CAOCS.  The  MM5  is  a 
limited-area,  non-hydrostatic,  terrain-following 
sigma-coordinate  model  designed  to  simulate  or 
predict  mesoscale  and  regional-scale 
atmospheric  circulation.  It  has  been  developed  at 
Penn  State  and  NCAR  as  a  community 
mesoscale  model  and  is  continuously  being 
improved  by  contributions  from  users  at  several 
universities  and  government  laboratories. 

2.2.2.  Land  Surface  Parameterization 

The  horizontal  resolution  is  60  km. 

Sixteen  pressure  levels  were  used  with  the  top  at 
10  mb.  We  use  a  soil  water  availability  function  to 
compute  the  soil  water  content  by  using  10 
specified  vegetation  types.  Dominant  deciduous 
and  coniferous  forests  cover  the  most  land  areas. 
There  are  small  areas  of  tropical  forest  in  Taiwan, 
and  Indo-China  Peninsula. 


2.2.3.  Initial  and  Lateral  Boundary  Conditions 

The  initial  and  horizontal  lateral  boundary 
conditions  for  wind,  temperature,  water  vapor, 
and  surface  pressure  are  interpolated  from 
analyses  of  observations  from  European  Center 
for  Medium-Range  Weather  Forecasts  (ECMWF). 
For  the  present  study,  they  were  projected  on  a 
spectral  T42  grid,  and  a  vertical  resolution  of  14 
pressure  levels.  The  initial  conditions  were  the 
fields  on  00Z  April  1,  1998.  The  lateral  boundary 
conditions  were  provided  via  a  relaxation  method 
and  updated  every  12  hours. 

2.3.  Ocean-Atmosphere  Coupling 

Atmospheric  and  oceanic  surface  fluxes 
of  water,  heat  (excluding  solar  radiation),  and 
momentum  are  of  opposite  sign  and  are  applied 
synchronously.  Flux  adjustments  are  not  used. 

We  use  the  same  flux  parameterization  as  the 
MM5  for  CAOCS.  The  only  difference  between 
the  two  is  the  use  of  sea  surface  temperature 
(SST)  data.  In  stand-alone  MM5  model,  SST  is 
prescribed  as  a  given  parameter.  However,  in  the 
CAOCS,  SST  is  predicted  by  POM  at  each  time 
step  for  the  atmospheric  component  (MM5). 

3.  PREDICTED  EAMS  CIRCULATION 

We  integrate  the  CAOCS  model  for  four 
months  from  1  April  1998.  Figures  1b,c,d  show 
the  predicted  temperature  and  velocity  vector 
fields  at  three  different  depths  (surface,  50m,  and 
200  m)  in  June  1998.  The  North  Equatorial 
Current  (NEC)  was  well  predicted  in  all  depths. 
After  encountering  the  western  boundary  along 
the  Philippine  coast,  the  NEC  bifurcates  into  the 
northward  flowing  Kuroshio  and  the  southward 
flowing  Mindanao  Current  (MC). 

The  MC  and  the  simulated  New  Guinea 
Coastal  Undercurrent  flow  equatorward,  feeding 
both  the  North  Equatorial  Counter  Current 
(NECC)  and  the  flow  from  the  Pacific  to  the  Indian 
Ocean,  that  is,  the  Indonesian  throughflow.  The 
MC  is  an  essential  part  of  the  shallow,  wind  and 
buoyancy  flux  driven  subtropical  circulation  cell 
that  carries  the  subtropical  waters  to  the  equator. 
It  is  typically  colder  and  fresher  than  the  South 
Equatorial  Current  (SEC)  and  the  New  Guinea 
Coastal  Current  over  the  thermocline,  due  to 
ventilation  in  the  North  Pacific  subtropical  gyre. 
The  northward  flowing  Kuroshio  is  the  major 
western  boundary  current  in  the  Pacific  Ocean. 
The  Kuroshio  Water  intrudes  into  the  SCS 
through  the  southern  Luzon  Strait  (the  strait 
between  Taiwn  and  Phillippines),  loops  around 


anticyclonically,  bifurcates  southwest  of  Taiwan 
into  the  northeast  flowing  South  China  Sea 
Warm  Current  (SCSWC)  and  southeast  flowing 
current  which  exits  the  SCS  through  the  northern 
Luzon  Strait.  In  the  SCS,  a  multi-eddy  structure 
was  predicted  with  central  SCS  warm  anticyclonic 
eddy  with  several  surrounding  cyclonic  eddies. 
The  continuation  of  the  SCSWC  becomes  the 
Taiwan  Warm  Current  after  passing  the  Taiwan 
Strait.  Major  current  system  in  the  EAMS  were 
realistically  predicted  using  the  CAOCS.  The  flow 
pattern  is  consistent  with  recent  observational 
studies  (e.g.,  Chu  and  Li  2000)  and  earlier 
studies.  Interested  readers  are  referred  to  a 
recent  review  paper  on  the  EAMS  (Su  1998). 

The  CAOCS  model  was  assessed  using 
the  data  collected  from  the  intensive 
observational  period  (IOP)  of  the  international 
South  China  Sea  Monsoon  Experiment 
(SCSMEX)  and  from  the  Navy's  AXBT 
measurements  in  April-August  1998. 

4.  CONCLUSIONS 

The  CAOCS  has  capability  to  predict  the 
current  system  and  thermal  structure  of  the 
EAMS.  Comparing  to  the  SCSMEX  and  AXBT 
measurements,  CAOCS  successfully  predicted 
the  South  China  Sea  multi-eddy  structure  during 
the  SCSMEX  IOP  in  1998.  These  warm-core  and 
cool-core  eddies  have  radii  varying  from  100  km 
to  300  km  and  maximum  tangential  velocities 
ranging  from  10  cm/s  to  20  cm/s.  The  cool-core 
eddies  are  cyclonic  and  the  warm-core 
eddies  are  anticyclonic. 
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